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Abstract 

An open quantum system consists of leads connected to a device of interest. Within the nonequi- 
librium Green's function technique, we examine the replacement of leads by self-energies in con- 
tinuum calculations. Our starting point is a formulation of the problem for continuum systems by 
T.E. Feuchtwang. In this approach there is considerable flexibility in the choice of unperturbed 
Green's functions. We examine the consequences of this freedom on the treatment of leads. For 
any choice the leads can be replaced by coupling self-energies which are simple functions of en- 
ergy. We find that the retarded self-energy Y7 depends on the details of the choice of unperturbed 
Green's function, and can take any value. However, the nonequilibrum self-energy or scattering 
function can be taken to be independent of this choice. Expressed in terms of these self-energies, 
nonequilibrium transport calculations take a particularly simple form. 

PACS numbers: 73.23.-b, 73.63.-b,73.40.-c 
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I. INTRODUCTION 



A general method for calculating electronic transport properties is provided by the use of 
non-equilibrium Green's functions (NEGF) .S'li'lll Because it has a fully quantum mechanical 
and non-equilibrium formulation, and because its foundational assumptions are transparent, 
this technique has been used to study a very wide range of transport problems. 

In this paper we will examine the way in which the leads permitting transport through 
an open quantum system (a "device") can be replaced by self-energy terms in a Dyson's 
equation. The analysis is limited to the case where the leads can be well-approximated as 
systems of non-interacting electrons in local thermodynamic equilibrium. In principle this 
approximation is not very restrictive, since sufficiently far even from a mesoscopic device 
the wires carrying current become macroscopic. Notice that the device whose transport 
properties are being calculated is left completely unspecified - it can contain interactions, 
disorder, etc. 

A similar treatment of the leads has been used by Datta, in an approximation in which 
the spatial continuum is replaced by a lattice of discrete points.! An apparently very dif- 
ferent approach for continuum systems was developed in a series of papers by Caroli and 
co-workers0 and by Feuchtwang.i The Feuchtwang approach has been the basis of NEGF 
calculations by, for example, Ref. y. Other workers have also cast the continuum problem 
into a self-energy form.ini 

Our point here is to show that the Feuchtwang approach for the continuum also can 
be repackaged so that leads are replaced exactly by simple self-energies, and to explore 
this possibility in some detail. We make two simple changes to Feuchtwang's calculation: 
regrouping the Dyson's equations to eliminate differential operators, and then using analytic 
continuation as provided by Langreth's theorem! (and learned by us from Haug and Jaug hJ; 
we will refer to this admirable textbook for all standard information about NEGF). The 
result is in our opinion pleasingly clear, and simplifies transport calculations of electronic 
devices. 

There is some freedom in how the self-energy is identified. A Hamiltonian H can in 
principle be arbitrarily separated into an unperturbed Hq and a perturbation H — Hq. If 
calculations can then be carried out exactly, the results will be independent of the initial 
separation. There is a similar arbitrariness in the definition of the unperturbed retarded 
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and its corresponding retarded self-energy S*". (In the Feuchtwang approach this comes 
from an arbitrary choice of internal boundary conditions). The exact Green's functions are 
of course independent of this choice. It turns out, moreover, that the "less-than" self-energy 
S*^ is independent of this choice, for a broad class of internal boundary conditions. 

For clarity we will focus on the simplest case: a one- dimensional system of non-interacting, 
spinless electrons in steady-state. But the results are rather general. For example, electron- 
electron and other interactions can be included by adding their appropriate self-energies to 
the coupling self-energy described here. 

Extension to two and three dimensions is also straightforward. One motivation for this 
paper is that the approach we describe turns out to work well even for systems with irreg- 
ular shapes.0 Because the final results are independent of the choice of internal boundary 
conditions, in practice one should make whatever choice is simplest. In one dimension any 
choice is equally easy, but in higher dimensions with irregular shapes, the right choice can 
simplify the calculation considerably.0 

The model we will consider is an open one-dimensional system consisting of two leads L 
and R (at X < L and x > R, respectively), connected to a device D (within L < x < R). 
This has Hamiltoniani 

H = e{L- x)Hl + e{x - L)e{R - x)Hd + e{x - r)Hr. (i) 

where 6' is a step function (0 to 1) and 

i^a = -|^^ + V;(x), a = L,D,R. (2) 

The single-particle potential energies Va{x) are arbitrary. To give definite examples in 
the following we will usually assume that Vl,r (but not Vd) are constants. The leads 
are treated as reservoirs of non-interacting electrons in local thermodynamic equilibrium 
described by local chemical potentials An applied bias eV would be included by 

setting = fJ'B. + and perhaps Vl = Vr + eV. This would lead to a current flow /, 
and a typical purpose of a NEGF calculation is to find the I-V characteristics of the device. 
The current and other physical observables of a non-equilibrium system are obtained from 
the total "less-than" function G^. We concentrate here on finding G^. The path we follow 
is to calculate the retarded Green's function G^ and then use analytic continuation rules to 
obtain G<. 
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II. RETARDED GREEN'S FUNCTION 



When the Hamiltonian is independent of time, the retarded Green's function C = 
0:2^2) can be Fourier transformed in ti — ^2- For brevity denote the resuh G{12) = 
G{xiX2Uj). For a system of noninteracting electrons this satisfies the equation of motion 

[huj - H{1) + ir/]G"(12) = 6{xi - X2). (3) 

Here H{1) is the total Hamiltonian in Eq. (P written in terms of coordinate xi, and rj > is 
infinitesimal. The advanced Green's functions, used below, are simply obtained by G"(12) = 
[^'■(12)]*. 



A. Unperturbed 

Following Feuchtwang,i one can similarly define individual "unperturbed" retarded 
Green's functions associated with each region a = L,D, R: 

[hu - H^{1) + ivK{12) = 6{x^ - x^), (4) 

subject to internal boundary conditions at a; = L, R. Like Feuchtwang we will consider the 
set of possible homogeneous boundary conditions 

^^£(12 



dxi 



= 7i^£(L2), ifx2< L (5a) 

xi=L 



dxi 
5^/£,(12) 



-fDL 90(^2), iiL<X2<R (5b) 



X^ =L 



dxi 



= jDRgniR'^), ifL<X2<R (5c) 
= 7ii^?fi(^2), ifx2>i?, (5d) 

xi=R 



dxi 

plus identical conditions with xi and X2 interchanged. Unlike Feuchtwang, we do not assume 
that the 7 are all equal. (This is important in higher dimensions.0) Eqs. (|^) include Dirichlet 
and von Neumann boundary conditions as the limiting cases 7 00, 0, respectively. The 
total Green's functions G cannot, and as we will see do not, depend on these internal 
boundary conditions. 

One also needs boundary conditions at the system's physical boundaries Xi — > ±00. These 
can be any homogeneous conditions. As long as we choose G to have the same homogeneous 
conditions on the physical boundaries, the boundary conditions there play no role. 



Both the perturbed and unperturbed retarded functions have matching conditions of the 
form 

,,12)1- _..0 (6, 

^■^^ xi=X2-e ''' 

for infinitesimal e > 0. The first of these comes from integrating Eq. across an infinites- 
imal interval about xi = X2- 

When Vl and Vr are constant, it is simple to obtain the unperturbed g^^ in the leads 
(a = L,R). Solve Eq. (|) separately in the regions Xi < X2 and xi > X2, and then use the 
boundary conditions Eq. (|^a,d) and the matching conditions Eq. (^) to fix the undetermined 
constants. The infinitesimal t] in Eq. (^) serves to pick out the outgoing (retarded) solution 
and can then be set to zero. The results are 

2m 



glil2) = J^e-^['=^(^<-^)+^^l sin[A;i(x> - L) + 0^] (7a) 
2m 



g'j,{12) = _ e^l^M->-K)+0«j sm[kn{x^ - R) + M (7b) 

where x> (a;<) is the greater (lesser) of Xi and X2, and where 

hu = V^ + (8a) 
2m 

7„ = kaCOt(f)a, (8b) 

for a = L,R. If ■jl^r are real, so too are (pn^i- It is apparent in Eqs. (0) that the unperturbed 
g^'{12) are symmetric under xi ^ X2- This is true for all other Green's functions as well. 



B. Total C 



Once more following Feuchtwang,! consider the following integral over lead L: 

dx3 {G''{32)[nuj - Hl{3) + ir]]gl{13) - gl{l3)[huj - H{3) + ir]]G'{32)} 

= e{L - xi)G^(12) - e{L - X2)gl{12). (9) 

The equality follows from the equations of motion, Eqs. (||,|). On the left-hand side of 
Eq. (P), all the terms in the two square brackets cancel except those involving derivatives. 
Hence the above integral also equals 



2m V ^ dxl 



9l- 



dxl 



2m 



dx^ 



dx-. 



(10) 



X3=L 



Here the right-hand side follows from Green's theorem. (There are no contributions from 
the physical boundary 0:3 —>■ —00 as long as (?£ and G have identical homogeneous conditions 
there.) Equating the right-hand sides of Eqs. (PJTO|), yields 



G%12) = gl{12)9{L-X2) + 



2m 



G'^(32)-^l(13) 



gG"-(32) 



{xi < L). (11a) 



We can perform similar integrals over the other two regions, replacing by (7^ or (yfjj, as 
appropriate, with the result: 



G'{12) = g}y{12)9{x2-L)9{R-X2) 



2m 



(9x3 



G^^m-gnm 



dG'X32) 



dx^ 



{L < xi < R), 



X3=L 



G^\12) = gUl2)e{x2 -R)-^ 



dx^ 



G'\32)-g^^iU) 



dG'{32) 

9X3 



(lib) 

(xi > /?j(llc) 



Each of Eqs. (11) is of the form of a Dyson equation (G"" = g^ + g^TI'G^'). They can be 
assembled into a single Dyson's equation for all space.l However, the self-energies include 
differential operators, an unattractive feature. 



C. Replace the leads by coupling self-energies 



Eliminating the derivatives gives a simpler expression. This is where we diverge from 
Feuchtwang's path. Let us focus on the only case of interest, G"'(12) when both arguments 
lie within the device (L < Xi,X2 < R). This is given by Eq. ([TT|b). Replacing dg]^/dx^ at 
the internal boundaries using Eq. (P3,c), 



G"-(12) = ^1,(12) + ^g'om 

- —gUlL) 



1drG'{R2) - 
1dlG'{L2) - 



9G"(32) 



9X3 

9G"(32) 



9X3 



X[i=R 



Xi=L 



(12) 



when L < Xi,X2 < R- Here we need dG^{32)/dx3 at X3 = L+. This is obtained from 
Eq. (|TT]a) by using Eq. (|a), choosing X2 > L, and letting xi ^ L^: 



9G''(12) 
9xi 



7l 



2m 1 



G"'(L2), X2 > L. 



(13) 
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The total C and its derivative are continuous across the internal boundaries. Consequently 
Eq. (|T3|) (and a similar expression for the derivative at R) give 



dG'{32) 



2m 1 



7L-— )G"'(^2), X2 > L (14a) 



2m 1 \ 



Plugging these derivatives into Eq. (|12D gives the following Dyson equation for in the 
device: 

G^(12) = gUl2)+g},{lL)J:i{u)G^{L2)+gUlR)^R{oo)G^{R2) 

(L < xi, X2 < -R), (15) 

where 

S2(-) = -|-(7.,-7.)-^, (16a) 

Thus for the purpose of finding G^ in the device, the contribution of the leads has been 
replaced by simple coupling self-energies ^. Notice that these are now c-numbers and 
not operators. This result is correct even for Dirichlet conditions as long as they are obtained 
as the limit 7^ — 00. One can get a final closed expression for G^{12) by evaluating Eq. (p!5[) 
at xi = L,R, solving the resulting pair of equations for G^' {L2) , G"^ {R2) , and substituting 
the results back into Eq. ([15|) . 

The standard case will have a constant potential Vl^r in the two leads, in which case we 
can substitute the particular expressions Eqs. (^ into Eqs. (p!6|): 

^liu;) = ^{-^j,^-rk^) (17a) 

^r{^) = iDR-ikR). (17b) 

We see that the choice of boundary condition in the leads drops out completely — there is 
no dependence on '^l,r left in Eqs. (|l^. However, the retarded self-energies Y7 do depend 
on the boundary conditions for the unperturbed g]^ in the device, via '-y^L and 715/2. We 
emphasize again the total G"" solved from Eq. (|1^) will be independent of this choice, also. 

These results agree with the discretized formulation by Dattai when the latter is taken 
to a continuum limit. Datta discretizes the coordinates into lattice points ja (with j an 



integer). In the Hamiltonian the kinetic energy term becomes a second difference. The step 
from the last site in lead L (at L — a, say) and the first site in the device (at L) becomes a 
perturbation. (A similar perturbation connects R in the device to i? + a in lead R.) Then 
an analysis like the above shows that the contribution of lead L to the discrete in the 
device can be written in terms of a self-energy 

TTj = -e'''^" (18) 



where 



huj = Al — coskLo). (19) 



ma? 



For small lattice constant a, ~ and Eq. (|T8|) becomes 



aEl^^{-a-'-ikL). (20) 



In the discrete approximation the unperturbed in the device vanishes at L — a. Thus for 
small a, 

(21) 



gjj{L - a, X2) = ^ gD\Lx2) - a 



xi=L 



dxi 

That is, the discrete calculation for small lattice constant a is equivalent to homogeneous 
boundary conditions with '-fci = a~^. Then the discrete self-energy Eq. ( pOl) agrees in the 
continuum limit with our general result Eq. (p!7| a). 

Because the boundary constants "ioL^lDR are arbitrary, one can obviously choose ^ 
with complete freedom. Probably the simplest case is to have von Neumann boundary 
conditions in the device (7_dl = = 7d_r). This choice is the only simple one in higher 
dimensions for devices of varying cross section.0 But every other choice is also possible. 
One can even choose /j to vanish (in which case 7/5^ = —ik^ and ''Jdr = ik^ are purely 
imaginary). The latter case amounts to building the contribution of the leads into the 
unperturbed (7^. 

III. NONEQUILIBRIUM GREEN'S FUNCTION G< 

In symbolic notation, the Dyson's equation Eq. ( plSj) for G"" in the device is 

G^- = g^^ + g'-j,T:G\ (22) 
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Standard analytic continuation rulesi then indicate that the exact contour-ordered Green's 
function G in the device satisfies 

G = gD + 9d^G, (23) 

where Qa is an unperturbed contour-ordered Green's function and S is the contour-ordered 
analogue of S^. Then a further application of the analytic continuation rules yields 

G^ = g^ + glT.^G< + g'j,i:<G'' + gp^'^G^ (24) 

Here G°-{12) = [G"'(12)]* is the advanced total Green's function and = (T,^)* the advanced 
self-energy. 

But what is the "less-than" self-energy S^? Our expressions Eqs. (p!6D for the retarded 
self-energies ^ involve inverses, and can't be analytically continued using the standard 
rules.i One could make the following physical argument J the self-energies considered here 
represent the coupling to leads which are in local equilibrium; hence one should be able to use 
the relationship between retarded and less-than functions valid for equilibrium quantities. 
That is, 

J:<{cu) = tUu;)T4u;) (25a) 

where 

r,(^) = ^K(a;)-S^(^)], (25b) 

for a = L,R. Here /a(ci;) = l/{exp[P{hLj — fia)] + 1) is the Fermi-Dirac equilibrium distri- 
bution for region a. As long as 'Jdl,1dr are real, substituting Eqs. (0) into Eqs. (^) gives 

j:<{u) = zU{u)^e{hu~V^), a = L,R. (26) 
m 

While this physical argument may not be entirely convincing, it has the merit of being 
correct (for real 7z)l,7z)_r), as we will now show. 

A. Unperturbed 

Before coupling the three regions, each is in local equilibrium with local chemical potential 
fia (a = L, D, R). Consequently we can use the equilibrium relationship between g^ and g^: 



g<{12) = zf4u;)A4l2) 
9 



(27a) 



where the spectral function in each region is: 



(27b) 



When all of the 7 in Eqs. are real, each satisfies boundary conditions exactly like 
those in Eqs. We will restrict ourselves to this case until the end of this section. 
When Vl^r are constant, the unperturbed g"^ in the leads (a = L, R) are 

5(<(12) = ifa{uj)9{huj - Va)-27^sm[ka{xi -a) + (pa] sm[ka{x2 -a) + (28) 

n rvrv 



using Eqs. ([7|, |27|) . The unperturbed g^ in the finite device is a sum over Dirac 5-functions. 
We will show later that it plays no role. 



B. Analytic continuation 



We calculate the total using the analytic continuation rules summarized byi 



(ABY 
{AB)< 



A'B' 

A''B< 



A^B" 



(29a) 
(29b) 



Here A (and B) are contour-ordered quantities, and are the retarded and less-than 



components of A. Begin with Eqs. (|Tl]). Using Eq. (|29| a), each of these can be analytically 
continued back to contour-ordered expressions. (Simply remove the superscript r). Then 
analytically continuing again with Eq. (p9[b) we have 



G<(12) 



' '5^?£(13)^<,,,, ^,,^,,9G<(32) 



+ 



8x3 

^^L(13)^a 



G<i32)-glil3)- 



8x3 

G<(12)^,<(12) + |- 



G^{32) - g<{13) 



8x3 

8G%32) 



8x3 



(xi < L,X2> L), (30a) 



8x3 



G<i32)-g},il3) 



X3=L 

8G< (32) 



8x3 



+ MMg"(32)-9S(13).^«°P2'^'' 



2m 



+ 



8x3 
8x3 
8x3 



G<(32)-<7le(13) 
G"^(32)-^7^(13) 



8x3 
8G< {32) 

8x3 
8G'' (32)' 
8x3 



{L<xi,X2<R), (30b) 



Xr,=L 



{xi> R,X2 < R). (30c) 



X3=R 
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The derivatives here can be ehminated much as they were in the Dyson equation for G^. 



Eqs. (pO|a,c) yield the derivatives dG'^/dx at L,R. The complex conjugate of Eqs. (|T3|) 
gives the derivatives dC"" / dx. Substituting all of these into Eqs. (pO|b) gives the following 
expression for in the device: 

G<(12) = g<{l2)+gl{lL)T.lG<{L2)+g^j,{lL)^<G\L2) 

+ ^^(1L)SJG"(L2) + {L^ R) {L < xi,X2 < R), (31) 

where SJ^^ are given by Eqs. (|1^) and their complex conjugate, and 

S<(cu) = ^"'"^ {a = L,R). (32) 

(?-(««)<(««) ^ ^ 



In the usual case the potential in the leads Vl^h is constant. Then, using Eqs. ([7| j27|) , this 
becomes equal to Eq. (p6[) . Thus the physical argument used at the beginning of this section 
is correct, at least as long as the 7's are real. 



Eq. (|3lD is of the standard form for the analytic continuation of a Dyson's equation for 
C to an equation for G^. It is easy to show that Eq. (|3TD can be written equivalently as 
the Keldysh quantum kinetic equationi 

G< = (1 + G"-S'-)^<(1 + S^G"^) + G"-S<G^ (33) 

For a finite length device, the first term in Eq. (^) vanishes. iJ^llli There are various ways 
to see this. It is clear on physical grounds that the exact Green's function including coupling 
to semi-infinite leads cannot depend on the initial occupancy of the finite device. We could 
choose fiD —* —00, in which case the device is initially unoccupied and the unperturbed 
g"^ vanishes. Proof that the term vanishes is based on the fact that the unperturbed spec- 
tral function for a finite device is a sum of Dirac 5-functions. The first term written 
symbolically in Eq. (Q) is, in detail, 

g<{l2) + G^{lL)T.lg<m + G^{lR)T:^g<{R2) 
+ g<{lL)i:iG^{L2)+g<{lR)i:iG^{R2) 
+ G^- {lL)^^^g<{LL)T,lG\L2) 
+ G\lL)i:ig<{LR)T.%G\R2) 
+ G^ {lR)Tr^g<{RR)T.lG\R2) 

+ G\lR)Tr^g<{RL)i:iG\L2). (34) 
11 



For the finite device we expand the unperturbed qds in discrete energy eigenstates ipii^x) of 
Hd with boundary conditions like Eq. (P3,c): 

= + , (35) 

i 

(7<(12) = -27ifnico)J2Mxi)i'iix2)5ihu - Si). (36) 

i 

Since each term in Eq. (|34|) contains a factor of g"^, the term can be nonzero only if hu equals 
some eigenenergy Si. But in this case [and solving for G"''"(L2), G^''^{R2) from Eq. (|15D], some 
algebra shows that the terms in Eq. (|3^) cancel exactly. 

Consequently the total NEGF is given by the second term indicated symbolically in 
Eq. (pSl). In detail for us this is 

G<(12) = G"(1L)S<G"^(L2) + G"(li?)E<G"(/22). (37) 

It seems from this last equation that the scattering functions ^ must be independent of the 
original choice of internal boundary conditions. After all, G"", G"", and are independent 
of this choice, and these three and are related by Eq. ([37|). However, this conclusion is 
not quite correct. We have shown that the are independent of the internal boundary 
conditions, and given by 

j:<{u;) = zUu;)^e{huj-V^), a = L,R, (38) 
m 

as long as the 7's are all real. However, when any of the 7 are complex, a similar analysis 
shows that pick up extra terms; but the extra terms give a cancelling contribution 
to Eq. (|37D . The accurate statement is rather that one can always choose the scattering 
functions to be of the form given by Eq. (^). 

IV. SUMMARY 

In this paper we have analyzed the replacement of leads by coupling self-energies in cal- 
culations of nonequilibrium Green's functions. This was based on a modification of Feucht- 
wang's approach for continuum systems. In general the form of a self-energy depends on the 
way in which the unperturbed Green's functions are defined. The final Green's functions 
of course do not. Here we investigated a class of homogeneous boundary conditions at the 
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junctions between device and leads. Because these internal boundary conditions are com- 
pletely arbitrary, the retarded self-energy E'' can take any value. The "less-than" self-energy 
or scattering function E^, on the other hand, can be taken to be entirely independent of the 
choice of boundary conditions. 

In one dimension no particular choice of boundary conditions is easier than another. 
In more complicated calculations, for example of two- and three-dimensional devices with 
nonuniform cross-sections, the ideas to take away from this work are: (1) The leads can be 
replaced by self-energies; (2) All internal boundary conditions give the same final result; so 
(3) Choose the boundary conditions which make the calculation easy. Following this route 
in our opinion simplifies transport calculations in general electronic devices of arbitrary 
cross-section. We will explore this in a subsequent paper.0 
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